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ABSTRACT 

We have performed a deep (35.5 hours exposure) spectroscopic survey of extremely 
red {I - K > 4) galaxies (ERGs) on the GOODS-South Field, using the GMOS 
spectrograph on the 8m Gemini South Telescope. We present here spectra and redshifts 
for 16 ERGs at 0.87 < z < 2.02, to a Hmit of Ks ~ 20.2. In the spectra of 10 of these 
galaxies we find emission lines of [OII]3727A, with fluxes corresponding to a mean 
star formation rate (uncorrected for dust) of 1.6 M0yr~^. For 3 ERGs we detect no 
emission lines and the remaining 3 lie at z > 1.5 where this line would be outside our 
wavelength range. Absorption features are seen in most ERG spectra. 

We examine the morphologies of these ERGs on HST-ACS images and fit radii and 
Sersic indices. We find three broad classes: (i) spheroidals; (ii) mergers at a variety of 
stages (some are spheroidal-disk pairs, some have tidal tails), often with a high surface 
brightness, and (iii) red spirals (which may have star-forming regions in their outer 
disks). 

We perform an age-dating analysis by fitting the spectra and 9-band photometry 
{BVIZJHK, plus 3.6/4.5 /im fluxes from Spitzer) of the ERGs with two-component 
models, consisting of passively evolving, old stellar populations combined with a 
younger, continuously star-forming component, for which the age and dust extinc- 
tion are allowed to vary from 10-800 Myr and E{B — V) = 0.0 to 0.5 mag. For only 
one ERG is the best-flt model purely passive, for the others the best-fit is obtained by 
including a star-forming component, which typically forms a few (0.26-13.5) percent 
of the stellar mass, and is subject to dust reddening averaging E{B — y) ~ 0.35. The 
ages of the star-forming components tend to be youngest (10-40 Myr) in the merging 
ERGs, and older (200-800 Myr) in spiral ERGs, with mixed ages for the spheroidals. 

The best-fitting mean ages for the old stellar populations range from 0.6 to 4.5 Gyr, 
averaging 2.1 Gyr, with masses from 3x10^*^ to 2xlO^^M0. The mean stellar formation 
redshifts of ERGs are spread from ~ 0.5 Gyr before the epoch of observation out to 
z ~ 5. We propose that most ERGs are galaxies, or mergers of galaxies, formed some 
Gyr earlier, in an early (z ~ 5) phase of massive galaxy formation, which since then 
have experienced a wide variety of merger and star- formation histories, accounting for 
the wide range of observed stellar ages. 

Finally, we examine the clustering of the ERGs on this field, as a function of the 
photometric redshifts estimated by Caputi et al. (2004). The comoving correlation 
radius is tq ~ 13/i^^ Mpc for the full sample and, dividing by redshift, is constant or 
increasing with redshift, thus favouring comoving (e = —1.2) evolution over a stable 
(e = 0) clustering model. 

Key words: galaxies: evolution, galaxies: high- redshift, galaxies: distances and red- 
shifts 
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1 INTRODUCTION 

The 'extremely red galaxies' (ERGs, also known as EROs) 
are a population of faint and distant galaxies, appearing 
in large numbers faintward of A" --^ 18, with very red 
optical/near- infrared colours {R — K > 5 or I — K > 4), 
corresponding to galaxies which are already old and pas- 
sively evolving at z > 1. The ERGs are of great importance 
in understanding galaxy evolution, in particular they repre- 
sent an evolutionary link between the rapidly star-forming 
galaxies detected by SCUBA at z > 2, and present-day mas- 
sive E/SO galaxies. 

A key property of ERGs is their strong clustering, as 
revealed by measurements of their angular correlation func- 
tion, ^(6*) (e.g. Daddi et al. 2000; Firth et al. 2002; Roche 
et al. 2002; Roche, Dunlop and Almaini 2003). Assuming 
a comoving model of clustering evolution, these observa- 
tions give for ERGs a correlation radius ro — 10-13/i~^ Mpc 
(where h — Ho/WO kms^^Mpc^^), which supports their in- 
terpretation as the progenitors of giant ellipticals. 

However, ERGs are actually a very diverse population 
of galaxies. Spectroscopic observations find that some are in- 
deed old, passive galaxies at z > 1 (e.g. Dunlop 1996); how- 
ever, at least half of the ERGs show emission lines (Cimatti 
et al. 2002; Yan et al. 2004a) indicating ongoing star for- 
mation. Radio, X-ray (Daddi et al. 2004) and mid-IR (Yan 
et al. 2004b) observations indicate that many emission-line 
ERGs are actually undergoing major starbursts, with star- 
formation rates (SFRs) as high as ~ 50-100 M0yr~^. Sec- 
ondly, deep X-ray surveys reveal at least ~ 10 per cent of 
ERGs to host obscured AGN (Alexander et al. 2002; Roche 
et al. 2003). Thirdly, ERGs are morphologically diverse; a 
mixture of bulges, disks, irregulars and mergers (e.g. Yan 
and Thompson 2003; Cimatti et al. 2003). 

Initially we (Roche et al. 2002) colour-selected a sam- 
ple of 158 ERGs with A" < 21 on the ELAIS N2 field. 
We found evidence of strong clustering and a mixture of 
morphologies, and identified 7 ERGs as radio sources (ei- 
ther AGN or powerful starbursts). We also fitted the ERG 
number counts with a model in which they evolve into the 
present-day E/SO population, through a combination of pas- 
sive luminosity evolution, merging, and a gradual increase 
with time in the comoving number density of red galaxies. 

For our second study (Roche, Dunlop and Almaini 2003, 
hereafter RDA03), deeper optical (HST-ACS) and near-IR 
(VLT-ISAAC) data enabled us to colour-select 198 ERGs to 
a fainter limit of 7^ = 22, on a central region of the Chan- 
dra/GOODS Deep Field South. For the clustering of these 
ERGs we estimated ro = 12.5 ± 1.2h~^ Mpc, and we iden- 
tified 17 ERGs with X-ray sources in the 1-Msec Chandra 
survey (Giacconi et al. 2002), from the X-ray fluxes diag- 
nosing 13 as obscured AGN and the other 4 as probable 
starbursts. 

In our third study, Caputi et al. (2004) fitted template 
spectra to the 7-band optical/near-IR photometry of these 
same ERGs to derive photometric redshift estimates, which 
in turn were used to estimate stellar masses and luminosity 
functions. Caputi et al. (2005) then extended this analysis to 
the full sample of if < 22 galaxies on this field (see Section 
6). 

In this paper we further investigate the CDFS ERGs by 
means of a deep optical spectroscopic survey, with the Gem- 



ini Multi-object Spectrograph (CMOS). We aim to mea- 
sure redshifts and emission lines for a sample of ERGs to a 
limit Ks = 20-21, and furthermore to estimate the ages of 
these galaxies by fitting evolving models to their spectra and 
broad-band fluxes (including new 3.6/4.5 pim fluxes from 
Spitzer). We shall interpret and discuss these results and 
make comparisons with the photometrically-based analysis 
of Caputi et al. (2004, 2005) and the recent spectrographic 
surveys of ERGs presented by Yan et al. (2004a) , McCarthy 
et al. (2004), Daddi et al. (2004 and 2005) Doherty et al. 
(2005) and Longhetti et al. (2005). 

In Section 2 of this paper we describe the optical/NIR 
data from RDA03 and our new GMOS observations, and in 
Section 3 the reduction and calibration of the spectroscopic 
data. In Section 4 we present spectra for the ERGs where we 
have secure or probable redshifts, fit models to estimate ages 
and masses, and measure emission-line fluxes. In Section 5 
we describe the morphology of the ERGs and estimate their 
size and surface brightness. In Section 6 we re-examine the 
clustering of the RDA03 ERGs as a function of spectroscopic 
and photometric redshift, and in Section 7 conclude with 
interpretation and discussion. 

Throughout this paper the assumed cosmology is flm ~ 
0.25, Oa = 0.75, and Ho = 70 km s~^Mpc~\ in accord with 
the recent (3 year WMAP) estimates of Spergel et al. (2006). 
This gives the age of the Universe as 14.16 Gyr. Magni- 
tudes are given in the Vega system except where labelled 
by subscript as AB magnitudes, the conversion between 

the two being {Ks,H, J,I775)aB = {Ks,H,J,l775)vega + 

(1.841,1.373,0.904,0.403). 



2 OBSERVATIONS 

2.1 Optical/NIR Imaging 

The Chandra Deep Field South forms one of the two fields 
of the publically available European Southern Observatories 
(ESQ) Great Observatories Origins Deep Survey (GOODS). 
As an early part of GOODS, deep imaging in the near-IR 
H , J and Ks bands was obtained for a central region of the 
CDFS, covering 50.4 arcmin^ and centered at RA 3''32™30^ 
Dec. -27 : 47 : 30 (equinox 2000.0), using the ISAAC cam- 
era on the 'Antu' Very Large Telescope (at the ESO site at 
Cerro Paranal, Chile). Also as part of GOODS, the central 
CDFS was mapped in four bands (-8435^606-^775-^850) with 
the Hubble Space Telescope Advanced Camera for Surveys 
(ACS) (Giavalico et al. 2004). Most recently, the GOODS 
Field South is currently being surveyed with Spitzer at 3.6- 
24 jj,m. For this paper we have available 3.6 and 4.5 fim 
fluxes for the galaxies on this field, derived from the 'epoch 
2' data release of May 2005. 

Firstly, to define our ERG sample, we ran the source 
detection program SExtractor on the mosaiced ISAAC Ks- 
band data, to obtain a sample of Ks < 22 galaxies (with 
completeness limit Ks — 21.5). For more details of this, see 
RDA03. With SExtractor we measured for these galaxies 
both total and aperture (2.0 arcsec diameter circular) mag- 
nitudes. 

These Ks < 22 galaxies were positionally matched to 
the corresponding detections on the ISAAC H and J mosaics 
and the ACS 1775 images. Colours were measured for each 
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Figure 1. Positions of tiic 198 ERGs in the CDFS sample of 
RDA03, within the borders of its 50.4 arcmin^ area. The Chan- 
dra X-ray sources (Giacconi et al. 2002) are also shown. Bold 
symbols indicate the 47 ERGs targeted by our GMOS spectro- 
scopic survey, and the 16 of these which form our redshift sample 
(see Section 3 and Tabic 1) are labelled with their ID numbers. 



as the differences in the 2.0 arcsec aperture magnitudes, and 
the ERG sample was then selected as the Ks < 22 galaxies 
redder than 7775 — Kg = 3.92, which numbered 198 in the 
50.4 arcmin'^ area. We numbered these objects 1 to 2555 
using their order of detection in the Ks < 22 catalog. 

For the galaxies studied in this paper we also measure 
3.6 and 4.5 /im fluxes from the Spitzer 'epoch 2' data release 
of May 2005. The Spitzer data has a lower resolution of 
FWHM ~ 1.6 arcsec. Because of this, larger apertures of 
2.8 arcsec diameter were employed to measure the Spitzer 
fluxes, and to these fluxes aperture corrections of 0.5 mag 
for 3.6/im and 0.55 mag for 4.5/im were applied, as derived 
from the point-spread function. 



2.2 Spectroscopic Observations 

The target galaxies for our GMOS survey form a subset 
of the RDA03 sample, consisting primarily of the brighter, 
Ks < 21, ERGs within a 5.5 x 5.5 arcmin area that can be 
covered by a single GMOS mask. A spectrograph mask was 
prepared with 47 slitlets (each of length 2.0 arcsec and width 
1.0 arcsec) devoted to our target ERGs. Figure 1 shows 
their positions within the larger area of our CDFS sample. 
Observations were performed with GMOS on the Gemini 
South telescope in Chile, over the period from 24 Septem- 
ber 2003 to 21 February 2004, using the low-dispersion, R150 
grating and central wavelengths in the region of 7500A. In 
total we were able to obtain 70 exposures, totalling 35.5 
hours, through the ERG slit mask. These science observa- 



tions were interspersed with CuAr arc lamp and standard 
star exposures, for wavelength and flux calibration. 

These observations utilized GMOS in the the 'nod-and- 
shuffle' mode described in detail by Abraham et al. (2004). 
The sky background is a serious problem for ground-based 
spectroscopy, especially at redder, A > 7500A, wavelengths 
where there are strong emission lines. The 'nod-and-shuffle' 
mode was designed to provide improved subtraction of the 
sky background, even where this is subject to considerable 
spatial and time variation. Essentially, the telescope is 'nod- 
ded' (by 1 arcsec in this case) along the spectroscopic slit, 
observing at each position alternately in a rapid one-minute 
cycle. In position A the object will be towards one side of 
the slit and the adjacent sky on the other, in position B the 
object will be on the other side. With each cycle, the charges 
stored on the CCD from the two positions are 'shuffled' (by 
2 arcsec) into adjacent strips on the CCD frame, where they 
are accumulated over the course of each half-hour exposure. 
In this way the sky is observed simultaneously with, and as 
close as possible in position to, each individual target galaxy. 



3 DATA REDUCTION 

GMOS data reduction was carried out using IRAF rou- 
tines, including some specifically written for this instrument. 
Firstly, bias frames were subtracted from the target galaxy, 
arc-lamp and standard star exposures. The next step was 
sky subtraction, performed individually for each of the 70 
exposures. The nod-and-shuffle mode produces a CCD frame 
showing two adjacent strips for each spectrograph slit, cor- 
responding to the two nod positions. To sky subtract, the 
frame is shifted by the shuffle step and subtracted it from 
itself (using the IRAF 'gnsskysub' task). This produces a pos- 
itive and a negative sky-subtracted spectrum for each target 
galaxy. 

Small dithers in the spectral and spatial directions be- 
tween the 70 exposures (introduced intentionally for the re- 
moval of charge traps; Abraham et al. 2004) were removed 
by registering the 70 exposures in both axes. The 70 sky- 
subtracted frames could then be combined, using IRAF 'im- 
combine' with a 'sigclip' rejection to remove most of the 
cosmic rays, into a single 35.5 hour image, from which we 
extract the spectra. 

IRAF tasks (e.g. 'apall') were then employed to extract 
the 2D spectra from each slit, summing each in the spa- 
tial direction to give a ID spectrum. With nod and shuffle 
it is necessary to extract both the positive and a negative 
spectrum for each object, and then subtract one from the 
other. 

The galaxy spectra were wavelength calibrated using 
the CuAr lamp calibration spectrum as observed through 
the same slit of the mask. Each pixel is approximately 3.5A 
and the resolution of our spectra FWHM ~ ISA Flux cal- 
ibration was performed using a co-added spectrum of the 
calibration standard star LTT1020, and the spectral energy 
distribution of this star as tabulated by Hamuy et al. (1994) , 
to derive (using IRAF 'sensfunc') a sensitivity function for 
the spectrograph (Fig 2a). 

Corrections were also applied to the spectra for the tel- 
luric (Earth's atmosphere) absorption bands, using our cali- 
bration spectrum of a white-dwarf standard star (these hav- 
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Figure 2. (a) Sensitivity of the spectrograph as a function of 
wavelength (as measured using calibration star LTT1020), and 
(b) sky noise within an extracted spectrum, calibrated in i<\. 



ing few intrinsic spectral features), EG131. This was divided 
by a fit to its own continuum, giving a function almost flat 
at unity except for four telluric absorption bands, centred 
at 6560, 6875, 7610 and 9075A. These were extracted, con- 
verted into magnitude units, and inserted into a function 
y{X) = 1 to give a telluric correction function. 

The ERG spectra were divided by this correction func- 
tion and by the sensitivity function, to give spectra cali- 
brated in Fv (in units of 10~^^ ergs cm~^s~^Hz~^, i.e. /iJy), 
and F\ (in units 10"^^ ergs cm~^s~^A^^), ready for further 
analysis. 



4 SPECTRA OF THE ERGS 

4.1 Defining the Redshift Sample 

On examining the GMOS spectra we found that our survey 
had a good rate of success for ERGs to Ks — 20, but failed 
to reach the fainter Ks = 21 limit of the target list. Of our 47 
spectra of ERG targets, 16 were galaxy spectra of sufficient 
quality that we were able to determine good or, at least, 
likely redshifts. 17 showed no significant signal, 13 had some 
signal in the continuum but we were unable to confidently 
identify any spectral features above the noise, and one was 
identified as a red Galactic star. The 16 'successful' spectra 
include 13 of the 18 target ERGs with Ks < 20, 3 of the 6 
target ERGs with 20 < Ks < 20.25 but none fainter than 
this. 

Redshifts were estimated by a combination of (i) fitting 
models to the continuum, by the method described below, 
(ii) identifying absorption features, and (iii) from emission 
lines such as [Oil] 3727A, where present. For galaxies with 
an emission line we determine both an emission (zem) and 



a continnum-fit/absorption-line (zc/) redshift. In all cases 
we found the two consistent within the ~ O.Olz statistical 
error on z^f. Where both are given Zem is likely to be more 
precise. 

Two of these redshifts, for the galaxies 1843 and 1158, 
are based on spectral features of poor signal-to-noise and 
must be regarded as less reliable than the other 14. However 
these two redshifts are supported by very good fits of our 
model spectra to the 7-band photometry. 

We therefore define a redshift sample comprised of these 
16 ERGs. The 16 redshifts range from z = 0.0868 to z = 
2.017, with 6 in a 'spike' at 1.180 < z < 1.225, which may be 
indicative of large-scale structure. One of the 16 ERGs, ID 
no. 1843, corresponds to a Chandra source. This is XID:253 
in the list of Giacconi et al. (2002), with a flux and hardness 
ratio (HR = 0.66 ± 0.11) indicative of an obscured AGN 
(sec RDA03). 



4.2 Comparison with Spectroscopy of Other 
Authors 

Some of the objects in our redshift sample were recently 
observed spectroscopically by other authors, in separate and 
independent studies, published during the writing of this 
paper - Vanzella et al. (2005), who used F0RS2 on the 
VLT, and Doherty et al. (2005), using DEIMOS on Keck. 
For objects in the overlap between these surveys, we make 
a brief comparison of the spectra. 

(i) Our object 1152 is in Vanzella et al. (2005) as GDS 
J033239.64-274709.1. The FORS2 spectrum gave z = 1.317, 
consistent with our z = 1.309, and similarly shows weak 
[Oil] emission, and absorption lines. 1152 is also object 2800 
in Doherty et al. (2005), observed with DEIMOS giving 
z = 1.318, with H&K absorption.. Neither spectrum goes 
far enough into the red to see well the Balmer lines. 
(ii) Object 1029 was observed as GDS J033235. 79-274734.7. 
The FORS2 spectrum gave z — 1.223, consistent with our 
z — 1.221, and is very similar to ours, with strong [Oil] emis- 
sion and some absorption features, including (weak) Balmer 
lines. 

(iii) Object 1412 is GDS J033235. 78-274627.5. The F0RS2 
spectrum gave z = 1.094, consistent with our z = 1.092, and 
like ours shows [Oil] emission. The 118 and H9 lines we find 
are not identified on the F0RS2 spectrum, but it extends 
further into the red to reveal H"f and H5 absorption. Thus 
on the basis of both spectra this seems to be a strong Balmer 
absorption galaxy. 

(iv) Object 1777 was not observed in the FORS2 survey, 
(v) Object 2536 is GDS J033231. 83-274356. 2. The F0RS2 
spectrum gave z = 1.550, slightly higher than our z — 1.513. 
Like our GMOS spectrum it features Mgll (2800A) absorp- 
tion, and extends a little further redwards, detecting an [Oil] 
emission line. 

(vi) Object 1064 is GDS J033230. 02-274726. 8. The F0RS2 
spectrum gave z — 1.553, a close match to our z = 1.548, 
and detects Mgll (2800A) absorption, and a weak [Oil] emis- 
sion line at the reddest extreme, 

(vii/viii/ix) Objects 1143 and 1565 and 2401 were not ob- 
served in the F0RS2 survey. 

(x) Object 1876 is GDS J033230.34-274523.6. The F0RS2 
spectrum gave z — 1.223, slightly higher than our z = 1.196. 
The FORS2 spectrum, like ours, shows Mgll (2800A) ab- 
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sorption, [Oil] emission and Balmer lines, extending furtlier 
redwards to Hj. 1876 is also number 2158 in the Doherty 
et al. (2005) DEIMOS sample, where it has z = 1.222 and 
H&K, H(5 absorption, very weak [Oil], and is classed as E + e 
(xi) Object 2246 is CDS J033229.82-274510.8. The FORS2 
spectrum gave z — 1.225, slightly greater than our z = 
1.196. Like our spectrum it shows [Oil] emission and H&K 
absorption, the spectra are very similar but ours has better 
signal-to-noise. 

(xii/xiii) Objects 1823 and 1955 were not observed in the 
F0RS2 survey. 

(xiv) Object 1762 is CDS J033227.02-274407.2. The F0RS2 
spectrum gave z — 1.127, slightly higher than our z = 1.102, 
and like ours shows moderate [Oil] emission and weak 
Balmer absorption. The authors mark H^ but it is barely 
visible. 

(xv) Objects 1843 and 1158 were not observed in the F0RS2 
survey. However, the AGN 1843 was observed with VLT 
FORS in an earlier survey of Chandra sources. Szokoly et 
al. (2004) initially proposed z = 0.484 on the basis of a sin- 
gle unidentified line. Zheng et al. (2004) claimed z = 1.89 
on the basis of a single line, this time supported by pho- 
tometric redshifts from COMBO-17 and hyper-2. However, 
for our 9-band photometry we obtain a hyper-z solution of 
z — 2.06, and our continuum is significantly better fitted 
for z — 2.017 than z = 1.89. Considering that photometric 
redshifts are uncertain by at least 8%, we retain z — 2.017 
as a provisional and approximate estimate. Near infra-red 
spectroscopy should detect multiple emission lines and give 
an unambiguous answer. 

In summary, out of our redshift sample of 16, 8 objects 
were also observed by Vanzella et al. (2005) and/or Doherty 
et al. (2005), who obtained very similar spectra similar spec- 
tra (although in two cases with a redshift discrepancy that 
might be larger than expected from statistical error). Seven 
of our redshifts appear to be new. Finally, for the z ~ 2 
AGN, more data is needed. 



4.3 Model-fitting 

We fit our ERG spectra, at 5800 < A < 9000A where the 
signal-to- noise is greatest, with a series of 40 template spec- 
tra from Jimenez et al. (2004). These represent passively 
evolving stellar populations, with a Salpeter IMF and so- 
lar metallicity, at ages from 0.001 to 14 Gyr. We exper- 
imented with models of lower and higher metallicity but 
found no indication that these fitted any better. Further- 
more the colours (Caputi et al. 2004) of our ERGs and the 
ultraviolet spectroscopy of de Mello et al. (2004) support the 
assumption that ERGs have approximately solar metallicity, 
even to z ~ 2. 

From the scatter between the 17 'blank' (i.e. no signal) 
slit spectra, we derived a function representing the noise per 
pixel produced by the sky background within the area of a 
single spectroscopic slit (Fig 2b). This is much worse at the 
red end due to strong sky emission lines and the declining 
detector response. 

We fit models to a combination of the CMOS spectra 
and the 9-band {BVIZJHK, 3.6/4.5/xm) photometry, us- 
ing a least-squares (minimum x^) method, with each pixel of 
the spectrum weighted using this noise function. The GMOS 



spectrum and the 9-band photometry make similar contri- 
butions to to the combined x^j which is minimized in the 
fitting. 

An initial series of fits as a function of redshift pro- 
vided the continuum-fit redshift Zcf for each ERG. We then 
fixed the redshift at this value and fit with a two-component 
model, this being the sum of a single-age passively evolv- 
ing model, of age Tpas, and a 'starburst' component of 
much younger stars, represented by the Jimenez et al. (2004) 
model, integrated with a constant SFR for a time interval Tsb 
prior to the time of observation. It is important to note that 
this component is not a single-age population, but contains 
stars of all ages from zero to Tsi,. The young component is 
subject to dust-reddening, following the extinction curve of 
Calzetti et al. (2000), parameterized in terms of E{B — V). 
The normalization of the young component is parameterized 
as fsb the fraction of the emitted flux it contributes in the 
rest-frame B-band, 4500A. Hence it is described by three 
parameters, T^b, E{B — V) and fst- 

For each ERG, we searched for the best-fitting model, 
by calculating a x^ for the GMOS spectrum, plus the 9-band 
photometry compared with a set of model spectra, forming 
a grid of 18000 points in the 4-dimensional {Tpas, fsb, Tsb, 
E(B — V)) parameter space. Tsb was allowed to vary from 
10 to 800 Myr and E{B - V) from to 0.5 mag {Av ~ 2 
mag). The best-fit model was taken as the position within 
this grid where x^ has the minimum value Xmin- 

Error intervals were estimated for the fitted parameters, 
as follows. Using our grid of models and their associated x^, 
we find the surface in parameter space, enclosing the best-fit 
model, on which x^ takes the value Xmin + 2.71(Xmin/-^d/), 
where Ndf is the number of degrees of freedom, approxi- 
mately 900 here. This surface defines an approximately 90% 
confidence interval. The error on a given parameter, e.g. 
Tpas, is then taken as the projection of this entire surface 
onto the Tpas axis {not just its intercept with the axis). This 
projected interval will then take into account any degenera- 
cies between the fitted parameters, which would produce an 
elongated surface and hence widen the projected error in- 
terval. 

The Jimenez et al. (2004) models are normalized in 
units of solar masses and so the normalizations of the model 
fits can provide estimates of the stellar mass in each galaxy, 
for both the passive {Mpas) and the young {Msb) stellar 
components. However, the model fits to the 9-band aperture 
photometry always gave higher normalizations than fits to 
the spectra (by factors 1.1-2.5), the reason being that a 2.0 
arcsec diameter aperture admits a greater fraction of the 
light from an extended galaxy than the 1.0 arcsec spectro- 
graph slit. We therefore give stellar mass estimates with the 
higher normalization given by fitting models to the aperture 
photometry. 

4.4 Spectra and Ages of the Redsliift Sample 

Figure 3 shows plots of our 16 GMOS spectra, the 9-band 
photometry and best-fit models. Spectra and models are nor- 
malized to fit the aperture (rather than slit) photometry. 

Table 1 lists co-ordinates of the 16 galaxies and Ta- 
ble 2 the best-fit parameters Tpas, Mpas, Msb and fsb, with 
approximately 90 per cent confidence interval error bars. 
Also given are a flux- weighted mean star- formation redshift. 
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Zmsf, which corresponds to the look-back time Tpas ear- 
lier than the redshift of observation, the masses of the pas- 
sive and starburst components in the best-fit model, plus 
Mfab{= j-j — ^ip — ) the fraction of the total stellar mass in 
the starburst/young component. Table 3 gives the 9-band 
aperture photometry, absolute magnitude Mb (AB system) 
in the rest-frame B-band, derived from the best-fit model 
normalized to the 2.0 arcsec aperture magnitudes, and Mbt 
an absolute magnitude with an approximate correction to 
total magnitude given by adding the (small) difference be- 
tween aperture and total K magnitudes. 

We find a wide variation in Tpas within the sample, from 
0.6 Gyr to 4.5 Gyr, with a mean of 2.1 ± 0.3 Gyr, and in the 
mass Mpas, from 3 x 10"'Mq to almost 2 x IO^^Mq, with 
a mean 1.1 x lO^^Af0. Note that we have performed two 
independent reductions of the GMOS data, and the other 
reduction was independently fit with both Jimenez et al. 
(2004) and Bruzual and Chariot (2003) models. These anal- 
yses gave age estimates closely similar to, and consistent 
with, those presented here. 

Only for one ERG was the best-fit obtained for a purely 
passive model, f^t = 0. For the 15 where x^ i^ smaller 
with fsb > 0, the best-fit star-forming component made 
up anything from 0.26 to 13.5 per cent of the total stellar 
mass (mean 3.5 ± 1.1 per cent) The best-fit ages and dust- 
reddenings for the star-forming components vary widely be- 
tween individual ERGs, and cover the full range allowed in 
this model, 10-800 Myr (mean 364 ±96 Myr) and 0-0.5 mag 
(mean 0.35 ±0.05). 

For 6 ERGs, the best-fit star-forming component is a 
short-term, < 100 Myr starburst, always with some dust 
reddening varying from E{B — V)— 0.2 to 0.5. For 9 others, 
the best-fit is with a longer period (100-800 Myr) of contin- 
uous star-formation, in 3 cases with little/no dust but in the 
other 6 with heavy extinction of E{B — F) ~ 0.5. In Section 
7.3 we discuss how age and recent SF history is related to 
morphology. 



4.5 Spectral Lines 

Ten of the 16 ERG spectra show a single emission line, which 
in all cases is assumed to be [OII]3727A. In most cases this 
identification is supported by the position of other spectral 
features. [OII]3727A emission is an indicator of ongoing or 
very recent (within ~ 10^ yr) star-formation, and according 
to Kennicutt (1998), a luminosity Loii ergs s^^ in this line 
corresponds to a SFR 1.4 x lO'^^Loii Moyr"^. 

For 3 ERGs we find no emission at 3727A, indicating 
that star-formation has ceased, is at a very low level, and/or 
is heavily obscured. One of these (2401) has a spectrum fit 
by a purely passive model, but the photometry of the other 
two (1701 and 1955) is best-fit with the inclusion of a heavily 
reddened young component, suggesting these are dusty post- 
starbursts. Three ERGs lie at z > 1.5 where [Oil] would not 
be visible in our wavelength range, but two of these have 
[Oil] visible at the red extreme of the F0RS2 spectra (which 
have a slightly redder limit), and all three have strong fiuxes 
at the blue end. 

Absorption features are prominent in almost all ERG 
spectra, the Mgll close doublet at 2800A and/or H±K (3934, 
3968A), and the ~ 4000A break is generally a prominent 



Table 1. RA and Dec co-ordinates of the 16 ERGs in our redshift 
sample; and total K magnitude (Vega system) 



ID no. 


R.A. 


Dec. 


Ktotal 




(deg.) 


(deg.) 




1152 


53.165199 


-27.785870 


19.14 


1029 


53.149200 


-27.792976 


20.21 


1412 


53.149109 


-27.774279 


19.23 


1777 


53.144367 


-27.735580 


18.87 


2536 


53.132652 


-27.732313 


19.60 


1064 


53.125080 


-27.790777 


18.98 


1143 


53.106606 


-27.787319 


20.23 


1565 


53.116936 


-27.768448 


19.77 


2401 


53.109440 


-27.764095 


18.93 


1876 


53.126465 


-27.756538 


18.24 


2246 


53.124275 


-27.752996 


19.68 


1823 


53.107304 


-27.741930 


18.61 


1762 


53.112560 


-27.735302 


18.66 


1955 


53.123837 


-27.736311 


18.91 


1843 


53.083546 


-27.746445 


20.00 


1158 


53.075893 


-27.785822 


19.78 



feature. indicating the presence of old stellar populations. 
Three emission-line ERGs (1412, 1823, 1876) show H8 and 
119 Balmer absorption lines indicating a large intermediate- 
age (10^-10^ yr) component. 

The [OII]3727A line fiuxes were measured from the 
GMOS spectra (using IRAF 'splot'), and multiplied by the 
ratio of the continuum flux in the 2.0 arcsec aperture to that 
in the GMOS slit, so as obtain a line flux more representa- 
tive of the entire galaxy. The [Oil] line fluxes are listed in 
Table 4, with errors calculated by summing the sky noise 
function (Figure 2b) over the resolution element (18 A) cen- 
tred on the observed line wavelength. From this flux and the 
spectroscopic redshift we derive the [OII]3727A luminosity 
of the galaxy, which is uncorrected for dust, and from this a 
star- formation rate SFRoii, using the relation of Kennicutt 
(1988). 

These SFRs, which again are uncorrected for dust ex- 
tinction, average 1.61 ± 0.37 Meyr^^. There are two ways 
in which we can estimate a very approximate dust-corrected 
SFR. Firstly, if the emission lines are assumed to be sub- 
ject to the same dust extinction as the continuum, then 
from the Calzetti (2000) reddening curve the [Oil] extinc- 
tion A3727 = 5.86E{B — V^). If we apply this correction to 
the SFRs in Table 4, using the value of E{B - V) best-fit 
to the spectrum of each galaxy, the mean SFR of the 10 
emission-line galaxies increases to 12.6 M0yr~^ 

Secondly, if we assume that the SFR and dust extinction 
in the star-forming component is truly constant, we can esti- 
mate SFRsb = Msb/Tsb- This will correct for dust, as this is 
taken into account in calculating Msb, but for an individual 
ERG this is very approximate due to the large error bars on 
Tsb- For the full sample of ERGs (but excluding the AGN), 
the mean SFRsb is 24 Moyr"^ 
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Figure 3. GMOS spectra and broad-band 

(ACS/ISAAC/Spitzcr) photometry of 16 ERGs, calibrated in F^- 
For each the upper plot shows the GMOS spectrum (solid line), 
our best-fit two-component model (dotted), and, where non-zero, 
the starburst component of the model. Labelled are the age of the 
model passive component, and the timescale T^f,, fractional flux 
contribution at \reat = 4500A, f^i, and dust reddening E{B — V) 
of the starburst component, and obvious spectral features. The 
lower plot shows a wider wavelength range with the 9 broad-band 
fluxes, the GMOS spectrum (solid) and the best-fit models as in 
the upper plot. Both plots are normalized to the 2.0 arcsec aper- 
ture fiuxcs. 
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Table 2. Rcdshift ^spec (from emission lines, 



where detected, otherwise from continuum fits, z^f); best-fit age of the passively 



evolving component (Tpas)', fraction /^i, of the flux at \restframe. = 4500A produced by the starburst component; age T^i, of the constant- 
SFR starburst; dust reddening E(B — V) of the starburst (magnitudes); mean stellar formation redshift (zms/) corresponding to the 



(normalized to the 2.0 arcsec aperture fluxes), in units of the solar mass Mq; Mf^f,, the starburst component as a fraction of the total 
stellar mass; Mg the total absolute magnitude in rest-frame B (AB system). Errors are approximately 90% confidence intervals. 



ID no. 

1152 

1029 

1412 

1777 

2536 

1064 

1143 

1565 

2401 

1876 

2246 

1823 

1955 

1762 

1843 

1158 



^spec 

1.310 
1.221 
1.092 
1.223 
1.512 
1.548 
1.254 
1.108 
1.200 
1.196 
1.195 
0.868 
1.180 
1.102 
2.017 
1.379 



-+0.1 

-0.6 

K).5 

-0.2 

■1+0.2 

''-0.1 



Tpas (Gyr) 

1.5^ 

1.0^ 

l.OJ 
2.0 ±0.2 
0.6 ±0.05 
1.0 ±0.15 

4.5J 

2.0l 

2.5^ 

1.5j 

2.0^ 

4.5" 



2.5 
1.5 
1.5_ 
3.5^ 



hi.i 
-2.1 

fO.2 

-0.6 

fO.2 

^-0.1 

r+0.3 

^-0.1 

-1.7 
-0.2 

-0.3 
0.1 
0.5 

+0.1 
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'-'■'^"-0.03 
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10+^5 
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o.ol°-i 
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-+o 

-0.1 




0.5 



0.5 



-0.05 



0.5^ 



-+o 

-0.05 

0.2 ±0.1 
n Q+0.05 

^•■^-0.1 

5+0 
'^•°-0.05 



^msf 

1 95+006 
J^-»'J_0.30 

1 68+°-2^ 

1 on+0.08 
^■•-"'-0.04 

2.08l°:f2 
1.78 ± 0.03 

.1+0.10 



2.04^ 



0.09 



+0 

0.1 
+0.2 

-0.0 





0.5 
0.0. 

'^•5-0.05 
5+° 

0.4 ±0.05 

0.4 ±0.05 

5+° 



7 90+00 
'■^■J_4.79 

1 85+°-l° 
9 07+0.17 
^■•^'-0.07 
-+0.16 
0.04 
+0.34 
0.11 
1.09 
1.25 
+0.16 
-0.21 
0.04 
0.19 
+0.16 
0.81 
+2.53 
-2.17 



1.75 
2.02 
2.98 
2.32 
I.6O: 
3.42 
4.59 



1.58 X 10^1 
3.23 X 10^0 
4.70 X 10^0 
1.08 X 10" 
6.17 X 10"' 
1.47 X 10" 
7.29 X 10"' 
3.87 X 10"' 
1.63 X 10" 
1.94 X 10" 
8.27 X 10"' 
1.91 X 10" 
1.43 X 10" 
9.35 X 10^0 
1.85 X 10" 
1.21 X 10" 



7.36 X 10** 
6.58 X 10** 
6.22 X 10° 
1.70 X 10^0 
8.20 X 10** 
1.31 X 10'^ 
1.89 X 10** 
3.38 X 10** 


5.10 X 10° 
2.64 X 10** 
2.36 X lO^o 
1.07 X 10^0 
8.89 X 10** 
3.63 X 109 
2.22 X 10° 



Mfsb 
0.0046 
0.0199 
0.1168 
0.1355 
0.0131 
0.0089 
0.0026 
0.0087 


0.0256 
0.0032 
0.1100 
0.0695 
0.0094 
0.0193 
0.0183 



* Host galaxy of obscured AGN - age/mass estimates may be affected by non-stellar component. 

Table 3. Nine-band (HST-ACS, ISAAC and Spitzer IRAC) photometry of the 16 ERGs in the redshift sample, with fluxes measured 
in fixed 2.0 arcsec diameter circular apertures (2.8 arcsec for Spitzer data) and expressed in the AB magnitude system, where rriAg = 
—48.60 — 2.5 log Fi, and l^Jy is m^s = 23.90. Absolute magnitude in rest-frame B (AB system), as estimated from the model fit to 
the aperture photometry (Mg), and with an approximate correction to total magnitude Mgj-- 



ID no. 


-B435(AB) 


Vu)<i(AB) 


/775(AS) 


^850 (AB) 


J(AB) 


H{AB) 


Ks{AB) 


3.6^m 


4.5/im 


Mb 


Met 


1152 


26.06 


25.13 


23.83 


22.88 


21.83 


21.43 


21.16 


20.21 


19.93 


-21.51 


-21.69 


1029 


25.85 


25.56 


24.58 


23.62 


22.79 


22.32 


22.11 


21.27 


21.00 


-20.43 


-20.49 


1412 


26.03 


24.94 


23.70 


22.84 


22.06 


21.65 


21.21 


20.29 


20.02 


-20.77 


-20.91 


1777 


{n.d). 


25.13 


24.07 


23.11 


22.06 


21.75 


20.90 


20.04 


19.79 


-20.95 


-21.14 


2536 


25.71 


25.04 


24.07 


23.47 


22.48 


22.10 


21.51 


20.79 


20.52 


-21.70 


-21.77 


1064 


25.14 


24.76 


23.92 


23.12 


21.78 


21.57 


20.94 


20.08 


19.81 


-22.06 


-22.18 


1143 


27.28 


25.47 


25.40 


24.78 


23.33 


22.76 


22.15 


21.37 


21.10 


-19.73 


-19.81 


1565 


27.08 


25.92 


24.50 


23.73 


23.11 


22.56 


21.69 


20.88 


20.61 


-19.79 


-19.87 


2401 


27.42 


25.71 


24.18 


22.97 


21.97 


21.23 


20.95 


20.43 


20.16 


-21.03 


-21.21 


1876 


25.64 


24.25 


23.20 


22.23 


21.29 


21.13 


20.39 


19.60 


19.31 


-21.75 


-22.06 


2246 


26.81 


25.81 


24.39 


23.41 


22.38 


22.22 


21.54 


20.90 


20.63 


-20.53 


-20.55 


1823 


27.07 


24.74 


23.21 


22.58 


21.78 


20.99 


20.53 


20.15 


19.88 


-20.86 


-20.94 


1955 


27.35 


25.69 


24.01 


22.81 


21.70 


21.52 


20.89 


20.06 


19.78 


-20.99 


-21.13 


1762 


25.21 


24.23 


23.29 


22.43 


21.87 


21.34 


20.74 


20.12 


19.86 


-21.06 


-21.06 


1843 


25.62 


25.16 


24.55 


24.15 


23.05 


22.36 


22.02 


21.27 


21.00 


-21.96 


-22.14 


1158 


{n.d). 


26.99 


25.09 


24.16 


23.46 


22.54 


21.85 


20.65 


20.38 


-20.35 


-20.58 



(n.d)= not detected (too faint) in this passband. 
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Table 4. Equivalent width of the [Oil] 3727A emission Hne, con- 
verted to restframe, the line flux (/on) and luminosity (Loil): 
and the star- formation rate derived from this (SFRqii = 1.4 x 
10~*^Lon M0yr~^), for ERGs in our redshift sample (these 
quantities normalized to the ERG flux in a 2.0 arcsec diameter 
aperture), la upper limits are given for galaxies where no [Oil] 
emission is seen. 



ID no. 


EW[OII] 


/on 


Lqu 


SFRoii 




A 


10-1* ergs 


10« 


Moyr-i 




restframe 


s-icm-2 


ergs s-i 




1152 


5 


8.40 ± 2.90 


9.21 ±3.18 


1.29 ±0.45 


1029 


46 


38.1 ±2.8 


35.1 ±2.5 


4.91 ±0.35 


1412 


9 


11.34 ±2.45 


7.88 ±1.70 


1.10 ±0.22 


1777 




<3.33 


<3.08 


< 0.43 


2536 


7 


7 


7 


7 


1064 


7 


7 


7 


7 


1143 


21 


7.38 ±2.13 


7.21 ±2.08 


1.01 ±0.29 


1565 


31 


15.85 ±2.34 


11.43 ±1.69 


1.60 ±0.24 


2401 




<2.41 


<2.12 


< 0.30 


1876 


5 


10.26 ±2.76 


8.94 ±2.35 


1.25 ±0.33 


2246 


12 


10.22 ±1.57 


8.91 ±1.37 


1.25 ±0.19 


1823 


10 


18.58 ±1.65 


7.26 ± 0.65 


1.02 ±0.08 


1955 




< 2.20 


< 1.86 


< 0.26 


1762 


8 


10.48 ±3.32 


7.45 ±2.35 


1.04 ±0.33 


1843 


7 


7 


7 


7 


1158 


40 


10.51 ±3.36 


13.13 ±4.20 


1.84 ±0.59 



?: Galaxy is at 2: > 1.5 where [OII]3727A would not be visible in 
GMOS A range. 

5 MORPHOLOGY AND SIZE OF ERGS IN 
THE REDSHIFT SAMPLE 

The HST-ACS images, in the stacked form used here (ver- 
sion VI. from the STScI), have a pixelsize 0.03 arcsec 
and resolution FWHM ~ 0.08 arcsec (~ 0.66 kpc at z — 
1.2). Figure 5 shows all 16 galaxies as they appear on the 
7775-band ACS image, and Figure 6 the B435-band images 
{^restframe ~ 2000A), which havc poorer signal-to-noise but 
provide extra information by highlighting regions of star- 
formation (SF). 

The appearance of each ERG is described briefly be- 
low (by 'high-SB' we mean an 7775-band central surface 
brightness of about 0.033 counts pixeP^s^^ — 21.7 AB mag 
arcsec"'^, a value typical of the spheroidal ERGs). 

1152 ([on] emitter) is a high-SB spheroidal surrounded 
by a much fainter ring (post- merger?). Only the central nu- 
cleus is visible in the B image. 

1029 ([Oil] emitter) appears in / to be a normal 
spheroidal, with a bulge profile. However, the B image shows 
spots within the galaxy which could be SF regions. 

1412 ([on] emitter) is a large (~ 1.8 arcsec diameter) 
disk with spiral arms and a small central nucleus. In B the 
nucleus cannot be seen but two brighter regions in the outer 
arms are just visible. 

1777 (no [on]) is a barred spiral of large diameter (~ 
2.4 arcsec) and relatively low SB. A nucleus and two spiral 
arms are clearly visible in I. It does not appear disturbed. 
In the B image the nucleus and arms cannot be seen and all 
that is visible is a tiny bright spot within the N spiral arm. 

2536 (2 > 1.5; UV-bright) is a merging pair of separa- 
tion ~ 0.25 arcsec. The NW component is smaller but higher 
in SB, the SE is irregular with bright regions on the W and 



SE sides. In the B band both galaxies are visible, but only 
the W and SE parts of the irregular component and the cen- 
tre of the NW galaxy are prominent. Both galaxies appear 
to contain bright SF regions within redder disks/spheroids. 
This ERG has the youngest Tpas of our sample combined 
with a young ongoing starburst. 

1064 (2 > 1.5; UV-bright) appears in the I band to 
be a high-SB elliptical. However, in the B-band the nucleus 
disappears and is replaced by bright regions on either side 
of it, forming a bar. This bar is very asymmetric and there 
is no sign of spiral arms, implying that this galaxy is not a 
barred spiral but a merger. 

1143 ([Oil] emitter) is a very small speroidal, without 
visible features (in both / and B) but appears to be a dusty 
starburst. 

1565 ([OH] emitter) is a complex system, with four sepa- 
rate intensity peaks. The southernmost galaxy has the high- 
est SB in I and could be a spheroid. Only one of the 4 com- 
ponents, an amorphous galaxy at the NW, is visible in B. 
This could be a multiple merger of a star-forming and 3 
passive galaxies. 

2401 (no [OH]) is a small, roundish, high-SB spheroid. It 
is not visible in the _B-band. The spectrum gives no evidence 
of ongoing or recent SF. 

1876 ([OH] emitter) appears to be a large, face-on spiral 
with a high-SB nucleus, and two arms reaching to a radius 
~ 0.9 arcsec. The galaxy is ringed by bright regions, presum- 
ably star-forming hotspots on the basis of their visibility in 
the B image, with gaps in the ring at the N and S. This 
galaxy has a strong colour gradient; the nucleus (r < 0.5 
arcsec) has a 'passive' colour of {B — I)ab = 2.49, but the 
ring (0.5 < r < 1.2 arcsec) has [B — I)ab = 1.39, a colour 
similar to the bluest galaxies in this sample (1029, 1064 and 
the AGN 1843) and indicating an observer-frame visible flux 
dominated by young stars. This very interesting object is a 
red spiral with ongoing SF, apparently of long duration, con- 
centrated in the outermost rim. 

2246 ([OH] emitter) is a small, roundish, high-SB 
spheroid. In the B band it looks even more smaller, and 
asymmetric, suggesting a central SF region. The SF appears 
to be prolonged and this is probably a post-merger. 

1823 ([Oil] emitter) is a merger of nuclear separation ~ 
0.3 arcsec, apparently of a spheroid (N) and a more irregular 
galaxy (S) with a lower J-band SB and three intensity peaks. 
In B the N spheroid is barely visible (passive?) but the S 
irregular is bright (starburst?). 

1955 (no [on]) is a high-SB spheroid with slight asym- 
metry in the outer regions, possibly due to perturbation by 
a smaller neighbouring galaxy separated by 0.8 arcsec sepa- 
ration. B morphology is similar but faint. 

1762 ([OH] emitter) resembles several of the disk merg- 
ers in the Arp catalog. The / image shows two nuclei sepa- 
rated by 0.34 arcsec, of which the larger (N) is at the centre 
of a partial ring or curved tail of radius ~ 0.4 arcsec, and the 
smaller (S) has a long (< 1.0 arcsec) straight tail. In the B 
band the system appears much more fragmented, the large 
nucleus is less prominent but the S nucleus is bright and 
parts of the tails are visible, together with a bright region 
in the SW. This is a merger with extensive SF, and some of 
the bright clumps in the tails might be new dwarf galaxies 
being formed in the interaction. 
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Figure 4. ACS /775-band images, from the GOODS data version vl.O, of the 16 ERGs in our rcdshift sample; (top row, left to right) 
1152, 1029, 1412, 1777, (2nd row) 2536, 1064, 1143, 1565, (3rd row) 2401, 1876, 2246, 1823, (bottom row) 1955, 1762, 1843, 1158. Each 
picture shows a 3.9 X 3.9 arcsec area, with an inverse greyscale with black corresponding to 22.27 AB mag arcsec"^. Galaxy 1158, due 
to its very low SB, is shown smoothed by a FWHM = 0.07 arcsec Gaussian and scaled up in intensity by a factor 2, and its ISAAC 
ftTs-band image is shown inset in the lower right (at a scale 0.4 that of its ACS image). 
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Figure 5. ACS iJ45o-band images, from the GOODS data version vl.O, of the 16 ERGs in our redshift sample; (top row, left to right) 
1152, 1029, 1412, 1777, (2nd row) 2536, 1064, 1143, 1565, (3rd row) 2401, 1876, 2246, 1823, (bottom row) 1955, 1762, 1843, 1158. Each 
picture shows a 3.9 X 3.9 arcsec area, with an inverse greyscale with black corresponding to 24.03 AB mag arcsec"^. Galaxy 1158, due 
to its very low SB, is shown smoothed by a FWHM = 0.07 arcsec Gaussian and scaled up in intensity by a factor 2. 
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Table 5. Half-light radius and Sersic index of the ERGs from 
model profiles fitted to HST-ACS Zg5o-band images. Galaxy 1876 
is fitted as a combination disk (n = 1) + bulge (n = 4), and 
galaxy/ AGN 1843 as a combination bulge + point-source (n = 0). 



ID no. 


Sersic Index n 


rui (arcsec) 


r-hl (kpc) 


1152 


4.8 


0.969 


8.43 


1029 


5.1 


0.330 


2.84 


1412 


1.2 


0.402 


3.39 


1777 


1.0 


0.543 


4.68 


2536 


1.4 


0.321 


2.83 


1064 


3.5 


0.360 


3.18 


1143 


2.4 


0.054 


0.47 


1565 


1.4 


0.384 


3.25 


2401 


5.8 


0.411 


3.53 


1876 


(85%) 1.0 


0.969 


8.32 




(15%) 4.0 


0.060 


0.52 


2246 


4.2 


0.261 


2.24 


1823 


1.2 


0.258 


2.05 


1955 


4.4 


0.807 


6.92 


1762 


1.4 


0.582 


4.92 


1843 


(60%) 4.0 


0.675 


5.91 




(40%) 0.0 


- 


- 


1158 


1.0 


0.864 


7.56 



1843 (> 1.5, Chandra source, AGN) has a central SB 
twice that of the other spheroidals in this sample. It ap- 
pears to be a bright point-source within a faint extended 
host galaxy visible in both / and B. 

1158 ([Oil] emitter) is of very low SB in the /-band. It 
is more prominent on the K-h&nA image. It appears to be 
an almost edge-on disk or ring some 1.6 arcsec in long-axis 
diameter, without a bright nucleus and without much asym- 
metry. In B it is even fainter and almost invisible, even with 
smoothing, but probably of similar size and shape. It is most 
likely a red spiral galaxy with the 'ring' appearence suggest- 
ing that most SF is in the outer regions rather than the 
nucleus (as in 1876?). Sizes (half-light radii, thi) were esti- 
mated for the ERGs by fitting the profile of each in the HST- 
ACS Zgso-band (approximately the rest- frame B-band) with 
the Sersic profile I{r) = lo ~ exp[{ — )~], using a method 
in which this was convolved with the detailed form of the 
HST-ACS point-spread function (see Floyd et al. 2004). In 
this equation n is the Sersic index, n = 1 for an exponential 
(pure disk) profile and n — 4 for a de Vaucouleurs bulge 
profile. The half-light radius r^i is a function of ro and n. 
Table 4 gives PSF-corrected rhi and n for the 16 ERGs; note 
that two galaxies were not well-represented by a single Sersic 
profile and have been fitted with two component models. 

Morphologically, these ERGs can be placed into three 
broad classes: (i) spheroidals of high Sersic index (n) — 
4.39 ± 0.40 (1152, 1029 1143, 2401, 2246, 1955, 1843), (ii) 
red disk or spiral galaxies, of low Sersic index (exponential 
profiles) (n) = 1.05 ± 0.05 and without much sign of distur- 
bance (1412, 1777, 1876, 1158), and (iii) merging or highly 
disturbed galaxies (2536, 1064, 1565, 1823, 1762) typically 
with intermediate Sersic indices (n) = 1.78 ± 0.43. 

Figure 6 shows the Sersic model r^i for the 16 ERGs 
plotted against the estimated total absolute magnitude Mb 
(the Mbt from Table 3). Also plotted for comparison are, 
firstly, a mean Vhi — Mb relation for local galaxies, as derived 
from the Driver et al. (2005) relation of B-band effective 
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Figure 6. Half-light radius from fitted Sersic model plotted 
against total absolute magnitude in Bj^b, for the 16 ERGs in 
the redshift sample (symbols indicate morphology). For the AGN 
1843 only the 60% extended component of the flux is included 
in the plotted absolute magnitude (giving Mb = —21.59). Also 
plotted are the mean size-luminosity relations for galaxies in the 
Millenium Galaxy Catalog, derived from Driver et al. (2005), and 
for early-type galaxies in the Sloan Digital Sky Survey (Shen et 
al. 2003). 



surface brightness to absolute magnitude for the Millenium 
Galaxy Catalog, converted to Hq = 70 and AB magnitudes 
{Bab = Bmgc — 0.13). This sample is composed predomi- 
nantly of spirals. Secondly, to represent local spheroidals we 
plot the best-fit size- luminosity relation of early-type (Sersic 
n > 2.5) galaxies in the Sloan Digital Sky Survey (Shen et 
al. 2003), with an approximate conversion Bab = r- + 1.17 
(Driver et al. 2005). 

Amongst spheroidal and disk ERGs, there is consid- 
erable scatter in the size-luminosity relation, but generally 
little, or moderate (< 1 mag) enhancement in SB relative 
to local galaxies of the same type. This would be consistent 
with the passive evolution of ellipticals which formed Gyr 
earlier, or with the exponetially declining SFRs of early- 
type spirals. One ERG, the disk 1158, has a significantly 
lower SB than a typical local galaxy and on the basis of 
its spectrum (Figure 3) this is likely due to heavy internal 
dust extinction, as well as a diffuse morphology. Merging 
ERGs, in contrast, show in their size-luminosity relation a 
marked (1-2 mag) SB brightening relative to local galaxies, 
especially in the case of the z > 1.5 mergers 2536 and 1064, 
presumably due to intense, short-term starbursts. 



ERG CLUSTERING VS. 
REDSHIFT 



PHOTOMETRIC 



We rexamine the clustering of ERGs in the CDFS sample 
using the photometric redshift estimates of Caputi et al. 
(2004, 2005), which will allow us to obtain a more reliable 
and less model-dependent estimate of their clustering prop- 
erties than given in RDA03. 

The angular correlation function, u!{d), is a measure of 
the clustering of galaxies as projected on the plane of the sky. 
It is related to the intrinsic clustering in three dimensions. 
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parameterised by the correlation radius ro, by means of an 
integration known as Limber's formula (see e.g. Efstathiou 
et al. 1991, RDA03) 

If the two-point correlation function, expressed in 
proper co-ordinates, ^(r) is represented by the simple model 

where 7 ~ 1.8 (observationally) and e represents the clus- 
tering evolution (e = is stable and e — —1.2 is comoving 
clustering), then Limber's formula gives uj{6) = A^O^^"'^^' , 
where 

" (l + z)''-(^+') 



A^ 



C^Vq 



-[{Nizy]dz/[ I Niz)dz] 

'0 



/o x^-^z)^ 

where x{z) is the proper distance and C^ = 3.679 for 7 = 
1.8. 

Hence in order to estimate ro from the observable Aui,sl 
model or observed N{z) is required, and for deep samples, 
the ro will also depend on e. However, we can also express the 
clustering measured at any redshift in terms of a comoving 
correlation radius rco, without any assumptions about e, as 
in comoving space 

S,{rc,z) = {rc/rcoy 

In this system comoving clustering would be a constant rco 
and stable clustering would be rco oc {1 + z)~3 . 

Previously (RDA03) we calculated the angular corre- 
lation function, (jj{9) of the CDFS ERGs.We obtained a 
~ 3a detection of clustering at the completeness limit of 
Ks < 21.5, which appeared consistent with the ERG 1^(6') 
measurements of Daddi ct al. (2000) , Firth et al. (2002) and 
Roche et al. (2002). By combining all these we estimated 
rco = 12.5 ± lAh^^ Mpc for ERGs. However, this estimate 
is dependent on a model N{z), taken from a 'merging and 
negative density evolution' model, fitted to the ERG number 
counts without any redshift data. Furthermore the cluster- 
ing evolution remains undetermined. 

We can now obtain a less model-dependent estimate of 
rco, and investigate e, by including redshift data. We assign 
to each ERG a redshift estimate, which for the 16 galaxies 
in our redshift sample is the GMOS spectroscopic redshift 
and for the remainder, the photometric redshift estimate 
from Caputi et al. (2004, 2005). Objects with Zphot > 4 
are assigned redshifts following Section 4.6 of Caputi et al. 
(2004), in which only two ERGs with Ks < 22 (one with 
Ks < 21.5) are accepted as probable z > 4 galaxies, and 
four of our original sample are reclassified as stars. 

Figure 7 shows the distribution of true/estimated red- 
shifts for the 175 ERGs brighter than the approximate com- 
pleteness limit Ks = 21.5. The photometric N{z) is quite 
close to the RDA03 model, especially at the peak, although 
more dispersed to lower and higher redshifts, presumably as 
real galaxies are more varied in their evolution. 

The angular correlation functions uj{6) were calculated 
for the CDFS ERGs using the same methods as in RDA03. 
Each oj{e) was fitted with the function 'A„(6I""*^ - 12.24)' 
(see RDA03), to obtain a power- law amplitude, A^i- The er- 
rors on the Au, were estimated using the same 'bootstrap' 
method as RDA03 (essentially repeating the analysis for a 
series of data catalogs with small regions excluded and find- 
ing the scatter in the resulting estimates). 

We first consider the full Ks < 21.5 sample of 



ERG X(z) for K<ai.5 
Photometric and spectroscopic 
Spectroscopic only 
Model trom HDA03 



s:i_ 



Redshift (photometric or spectroscopic) 



Figure 7. The distribution N{z) of photometrically estimated 
(Caputi et al. 2004, 2005), or spectroscopie (this paper) redshifts 
for the 175 Ks < 21.5 ERGs in our CDFS sample (solid his- 
togram); spectroscopic redshifts only (dashed histogram) and the 
best-fit model ERG N{z) from RDA03 (dotted line). 



175 ERGs. From its u}{6) amplitude measurement (Table 
6) and photometric/spectroscopic N{z) (Figure 7), Lim- 
ber's formula gives the comoving correlation radius rco ~ 
13.13+2:i4h"^ Mpc. This is consistent with the RDA03 es- 
timate based on a model N{z). Alternatively if we set e = 
this result can be fit with a stable clustering model with 
ro = 23.9ll4:t4h-i Mpc (at z = 0). 

We then recalculate u}{6) for the ERGs divided into two 
subsamples in phot/spec redshift. For both the low and high 
redshift subsamples, the u:>{0) and N{z) are used to derive 
a comoving rco (Table 6). Dividing the sample reduces the 
significance of clustering detection, but in most cases it re- 
mained > 2(7. We repeat this analysis for redshift divides at 
z = 1.4, 1.6 and 1.8. 

The subset uj{6) amplitudes are given in Table 6 and 
plotted on Figure 8. It can immediately be seen that the 
mean of the rco we derive for the low and high-redshift sub- 
samples is a little lower than that of the full A"-limited sam- 
ple, by about 19 per cent, which corresponds to 32 per cent 
in the A^i- The reason for this is that photometric redshifts 
are not always accurate and thus the low and high redshift 
subsamples will both experience some contamination from 
galaxies in the other redshift range, which are not correlated 
with the correctly selected galaxies and dilute the observable 
clustering. If, for example, the true io{9) of the lower redshift 
ERGs is At and of the high Ah, and we examine a low Zphot 
subsample with a contamination a, the io{9) we calculate 
will be Aobs — (1 — a)^Ai + oP Ah- The observed shortfall 
could be accounted for by a ~ 0.20. We have not attempted 
to correct for this is deriving a rco, and thus the low and high 
redshift rco will both be slight underestimates, (by similar 
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Table 6. Angular correlation function {u){9)) amplitudes (given 
at 6 = 1 dcg. for a fitted d~'^-^ power-law), estimated for CDFS 
ERGs, to a limit Ks = 21.5 (Vega system), divided into two 
subsamples in photometric/spectroscopic redshift. This is shown 
for redshift divides 2 = 1.4, 1.6andl.8 



^phot/ spec 

range 


^mean 


No. of 
ERGs 


u){9) amp. 
10-3(deg0«) 


rcO 
fe-i Mpc 


All 

2 < 1.4 


1.742 
1.125 


175 
72 


6.54 ±2.02 
10.94 ±3.89 


in 01+2.12 

id.Ji_244 


2 > 1.4 
2 < 1.6 
2 > 1.6 


2.187 
1.231 

2.458 


103 
101 

74 


7.39 ±4.22 
6.61 ±2.93 
8.34 ± 3.46 


11.39+^f, 
12.321^:?? 


2 < 1.8 


1.273 


111 


6.42 ± 2.89 


q nS+203 
^■^•^-2.56 


2 > 1.8 


2.577 


64 


12.41 ±4.70 


i4n4+2-'^4 

i4.U4_3 2g 
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Figure 8. Comoving correlation radius rcO, derived from the uj{0) 
amplitude and photometric A'^(2),for a full Ks < 21.5 sample of 
ERGs, and for these ERGs divided by photo-2 into low and high 
redshift subsamples (dividing at 2 = 1.4, 1.6 and 1.8) plotted 
against sample mean redshift. Also plotted are comoving (e = —2) 
and stable-clustering (e = 0) models normalized to the full sam- 
ple's clustering. 



amounts, as here the two subsamples have similar size and 
clustering) . 

Secondly, it can be seem that the higher redshift sub- 
samples have a comoving ro similar to or even a little greater 
than the lower redshift ERGs. This is clearly more consis- 
tent with comoving clustering rather than stable clustering. 
The slight increase in r^o with redshift, if real could reflect 
a stronger clustering for more massive/luminous ERGs. 
Comoving clustering of ro ~ 12-13/i^^ Mpc is consistent 
with the clustering of both shallower (Daddi et al. 2000) 
and deeper ERG samples, e.g. Daddi et al (2004) estimate 
ro = 9-17h~^ Mpc at 1.7 < 2 < 2.3. Extrapolated to lower 
redshifts, this model would overpredict the clustering of local 
passive galaxies of all luminosities, ro — 7.2h^^ Mpc (Madg- 
wick et al. 2003) , but be consistent with the clustering of the 
most luminous {Mr < —22.27) early-types at 0.3 < z < 0.9, 
ro = 11.2± l.O/i"^ Mpc comoving (Brown et al. 2003). This 
can be explained if (i) many of the z > 1 ERGs undergo 
further merging to become very massive ellipticals, (ii) only 
some (~ half) of the moderate-luminosity E/SOs found lo- 



cally have evolved directly from z > 1 ERGs and the others 
are formed at 2 < 1 from mergers of less strongly clustered 
spirals, thus diluting the clustering. 

To better understand the evolution of ERGs it will be 
useful to compare their clustering with that of their their 
likely progenitors at higher redshifts, such as sub-mm galax- 
ies, for which the clustering will be measured in the ongoing 
SHADES survey. There is already evidence that the most 
massive Lyman break galaxies at z ~ 4 have a similar co- 
moving ro (11.4 ± 2h'^ Mpc) to the z ~ 1-2 ERGs (Allen 
et al. 2005). 



7 DISCUSSION 

7.1 Ages of ERGs 

We perform spectroscopic age-dating of the ERGs by fltting 
their SEDs with passively evolving models (from Jimenez 
et al. 2004), and obtain mean stellar ages ranging widely 
from 0.6 to 4.5 Gyr, with a mean of 2.1 ± 0.3 Gyr. We can 
compare these results with similar age-dating analyses per- 
formed by other authors, though must be careful to con- 
sider the differences in how the sample are selected. Firstly 
Daddi et al. (2004) give ages (from FORS2 spectra) for 9 
A'-selected {Ks < 20) very massive star-forming galaxies at 
1.7 < 2 < 2.3 - note that not all of these are red enough 
to be called ERGs. Daddi et al. (2005) give ages from low- 
resolution ACS spectra of 7 ERGs at 1.39 < 2 < 2.47, of 
which 5 are passive ellipticals and 2 weakly star-forming. 
McCarthy et al. (2004) using Gemini Deep Deep Survey 
spectra, give ages for a varied sample of 20 red {I — K > 3.5) 
galaxies at 1.3 < z < 2.2. Finally, Longhetti et al. (2005) 
using low-resolution NIR spectra (observed with NICS on 
the La Palma TNG), estimate mass- weighted ages for 10 
of the most bright {K' = 16.6-18.4) and massive ERGs 
{R- K' > 5.0) present at 2 ~ 1.5. 

Figure 9 plots stellar age Tpas against spectroscopic red- 
shift for all these samples, with loci of the redshifts corre- 
sponding to time Tpas before observation. 

For the first four samples the mean stellar ages of ERGs 
at all 1 < 2 < 2.5 are spread over the range between 2^5 
and shortly (a few x 10* yr) before the epoch of observation. 
They seem to divide into two groups with z^sf ~ 2 and 

Zmsf ~ 3-4. 

The Longhetti et al. (2005) sample of very luminous 
ERGs differs from the others in that it seems to be con- 
centrated towards higher formation redshifts - it contains 
5 galaxies with Zmsf > 5. These have stellar masses 3.8- 
9.4 X 10^^ Mq, more than twice the mass of any ERG in 
our sample. Hence their ages may be evidence that the very 
most massive spheroidals formed significantly earlier and/or 
formed their stars most quickly, as in the model of Granato 
et al. (2004). 

It is also notable that the formation redshifts appear to 
trace the plotted 2ms/ loci. The ages of the oldest galaxies at 
each redshift, as a function of redshift, is an important tracer 
of the expansion of the universe and hence the cosmological 
model (e.g. Jimenez and Loeb 2002). We would expect Zmsf 
for the oldest galaxies to be constant with observed redshift, 
and the agreement with the plotted loci would support a 
A ~ 0.75 cosmological constant (w = — 1) model. 
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Figure 9. Mean stellar age Tpas, derived from model-fitting to 
spectra, against redsliift for our ERG sample, the ii'-selectcd 
galaxies of Daddi et al. (2004), and the ERGs of Daddi et al. 
(2005), McCarthy et al. (2004) and Longhetti et al. (2005). Plot- 
ted loci show the redshifts corresponding to time Tpas before the 
redshift of observation. 



As Tpas represents a fiux-weighted mean stellar age, 
the true age of a galaxy since formation must always be 
greater, to a degree dependent on the detailed form of the 
star-formation history. One interpretation may therefore be 
that most or all of the red galaxies observed as ERGs at 
z = 1-2 are galaxies, or mergers of galaxies, which formed 
at 2 > 4-5, and that the spread in their stellar ages up to 
this maximum redshift can be attributed to differences in 
their individual evolution, i.e. the mergers and starbursts 
they have experienced. 

Caputi et al. (2004, 2005), using photometric redshift 
estimates derived from seven-band photometry of ERGs and 
other 7^ < 22 galaxies in this field, found the comoving num- 
ber density of massive (> 5 x lO^^M©) galaxies in a full K- 
selected sample to fall only slowly from z = 1.75 to z = 3.5, 
where a significant fraction, 20-25 per cent, of today's mas- 
sive galaxies were already in place. However, the comoving 
number density of the reddest subset of these galaxies, the 
ERGs, evolves much more rapidly, steadily increasing with 
time from z = 3.5 to z = 1, and following a trend which 
extrapolates to the present-day E/SO population at 2 = 0. 
At 2 = 1-1.5, the rising comoving number density of ERGs 
approached that of all if-selected galaxies at the higher red- 
shift of 2 ~ 3.5. 

Hence, this suggests that a substantial population of 
massive galaxies formed at 2 > 4 and were initially star- 
bursting and relatively blue, but with increasing age many 
became red enough to be classed as ERGs. Interestingly, 
for the brightest (/ < 24.5) Lyman break galaxies at 
2 ~ 4, Allen et al. (2005) measure strong clustering of 
ro = 11.4±2h^QQ Mpc comoving, consistent with these evolv- 
ing into the ERGs ai z — 1-2. 

To produce the evolution in ERG number density these 
galaxies must have entered the ERG class over a wide red- 
shift range 1 < z < 4, again implying much variation in their 
star-formation histories. The formation of so many massive 
galaxies at z > 4 places significant constraints on models of 
galaxy formation (e.g. Fontana et al. 2004) and may favour 



models with feedback between the star-formation in forma- 
tive spheroidals and the growth of supermassive black holes 
within them (Granato et al. 2004; Silva et al. 2005). 



7.2 Ongoing Star-formation Activity in ERGs 

Yan et al. (2004a) and Doherty et al. (2005), found that 
> 75 per cent of ERG spectra have prominent absorption 
features and 4000A breaks, while > 50 per cent also show 
some [OH] emission. We find the same for our slightly deeper 
sample. For only 3 ERGs do we find neither [OH] emission 
nor any indication of a young stellar population from a blue 
excess in the SED. This is consistent with the Doherty et al. 
(2005) estimate that 28 per cent of ERGs show no evidence 
of recent star-formation. 

For a clear majority of ERGs, 10/13 at 2 < 1.5, we 
detect a [OH]3727A emission line. This includes two galaxies 
where the line only just detected with an equivalent width 
~ 5A, while at the other extreme we have 3 ERGs with [OH] 
equivalent widths > 30A. For the 10 emission-line galaxies, 
the [OH]3727A fluxes - uncorrected for dust - correspond 
to a mean SFR of 1.63 MQyr"^, intermediate between the 
2.8 Moyr"^ mean SFRoii of Yan et al. (2004a) and the 1.1 
Moyr"^ of the Doherty et al. (2005) ERGs. 

Of course the true SFRs of the emission-line ERGs will 
be higher. On the basis of the dust extinction estimated 
from our model fits, we estimate the real mean SFR will be 
about 12-24 Moyr^^. This may still be an underestimate 
as emission line regions of current SF may suffer more ex- 
tinction than the continuum produced by SF over a longer 
period. Spitzer mid-IR observations (Yan et al. 2004b) of 
ERGs similarly selected to those in our sample indicate that 
about half have SFR > 12M0yr~^, and for these the mean 
SFR ~ 6OM0yr^^. We should be able to obtain genuinely 
dust-independent measurements of the SFRs for the indi- 
vidual ERGs in the GOODS-S from the deep Spitzer 24/im 
survey of this field scheduled for later this year. 



7.3 Morphology and Star-formation History 

Examining the HST-ACS images of ERGs we find ERGs to 
be a mixture of spheroidals, spirals and merging systems. 
The classifications by eye are confirmed by Sersic model fits 
to the radial profiles (Table 5). We discuss the nature and 
SF activity of the ERGs within each of these morphological 
categories. 



7. 3. 1 Spheroidals 

We class 7 ERGs, almost half the sample, as spheroidals. 
Spectroscopically they are varied, and this may be related to 
whether or not they have recently undergone a merger event 
or the accretion of a gas-rich galaxy. One spheroidal (2401) 
appears to be old and purely passive, with no evidence of SF 
in the last 2.5 Gyr. Another (1955) shows no [OH] emission 
but its near-IR colours suggest it may be a dust-reddened 
post-starburst galaxy. One of the spheroidal ERGs hosts an 
obscured AGN, it may also be star- forming. One is old with a 
young, dusty starburst (1143). The final three (2246, 1152, 
1029) are relatively young (1-2 Gyr) with the model fits 
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suggesting prolonged, low-level, non-dusty SF; and appear 
to show post-merging features. 

7.3.2 Mergers 

Of these 16 ERGs 5 are obvious mergers, with multiple nu- 
clei, and all of these are star- forming galaxies, on the basis of 
[Oil] emission. The two z > 1.5 mergers also have a high sur- 
face brightness (Figure 6). Presumably the SFRs and dust 
content are both increased by the interactions. Some ap- 
pear to be mixed mergers of passive spheroidals with bluer 
irregulars (e.g. 1565 and 1823), and one (1762) has promi- 
nent tidal tails. On the basis of our model fits, four of the 5 
mergers are relatively young (0.6-2.0 Gyr) and undergoing 
short-term (10-40 Myr) dusty starbursts. Notably, in our 
sample; 4 of the 5 galaxies with Taf, < 50 Myr are merg- 
ers (the other is an AGN). This would be expected from 
the current interaction acting to tidally trigger an intense 
but short-term starburst. The fifth merger, a close irregu- 
lar/spheroidal pair, (1823) difi'ers in that it is ~ 4.5 Gyr old 
with a long T^i, and Balmer absorption. This may be ex- 
plained if it is a late-stage merger where most starbursting 
occurred at earlier encounters of the two galaxies, and the 
elliptical is already very old. 

7.5.5 Spirals 

Perhaps a more surprising discovery is that of a third mor- 
phological class of ERGs which are neither ellipticals nor 
mergers, but red spirals. We find four of these. One ex- 
planation for the red colours of these galaxies is that they 
are very dusty, the ERG colour criterion selecting out the 
spirals at the very top end of the range of internal dust 
extinctions, the other is that they long ago ceased form- 
ing stars. According to our model, the four spirals differ 
widely in age (1.0-3.5 Gyr) but in all cases contain a young 
component with a long Tgj, (200-800 Myr) and heavy dust 
extinction {E{B — \^) ~ 0.5). Three show [Oil] emission, 
the fourth (1777) may therefore have ceased forming stars. 
Strong Balmer lines are seen in two of the spirals, 1412 and 
1876, and this reinforces the evidence from model fits that 
spiral ERGs tend to contain large intermediate-age (0.1-1.0 
Gyr) populations. 

It is also notable that spiral 1876 has a luminous red 
(passive?) nucleus ringed by blue star-forming regions, pro- 
ducing a strong colour gradient. Indeed, for 3 of the spiral 
ERGs, a comparison of the / and B images shows that (in 
contrast to the merging ERGs), their star-formation is not 
predominantly nuclear but in the outer disks. The z — 1.55 
ERG #4950 of Daddi et al. (2005), classed as type Sa, also 
has a red nucleus with a blue outer ring. 

The existence and spectral properties of very red spirals 
can be explained if they are galaxies with continuous star- 
formation which has been gradually slowed or completely 
truncated. This explains the long Tsb, Balmer absorption, 
and the absence of emission lines in one spiral (already trun- 
cated?). In addition to this, they are very dusty, as the ERG 
colour criterion would select out the spirals at the very top 
end of the wide range of internal dust extinctions. 

Very red spiral galaxies also occur at lower redshifts, 
where they make up some 0.3 per cent of all galaxies (Ya- 
mauchi and Goto 2004), and tend to be found in the outer 



regions of clusters (Couch et al. 2002; Goto et al. 2003). 
In the model of Bekki, Couch and Shioya (2002) a spiral 
galaxy, infalling into a cluster, experiences a truncation of 
star-formation activity, as its gas content is removed through 
ram-pressure stripping by the intra-cluster medium and/or 
by cluster tidal fields. Following this the galaxy would be- 
come very red in ~ 1 Gyr but could retain a spiral ap- 
pearence, with visible arms, for > 3 Gyr, although would ul- 
timately evolve into an SO type. This mechanism could give 
rise to red spirals at ERG redshifts, as both spiral galaxies 
(Dawson et al. 2003) and large ERG-rich clusters (e.g. Kurk 
et al. 2004) are already in place at z > 2. However, some spi- 
ral ERGs e.g. 1876, must be in relatively early stages of this 
process as much SF is still ongoing in their outer regions. 

More information on the nature of spiral ERGs and the 
relative importance of dust and 'truncation' will soon be 
available from Spitzer 24/im data, giving a dust-independent 
SFR, and spectroscopy at l-2jj,m, showing Balmer lines and 
Ha. It is likely that some sort of environmental effect is 
involved because ERGs in general are very much more clus- 
tered than typical spirals. 



7.4 Further investigations 

Ongoing surveys of the CDFS with the Spitzer infra- 
red telescope will, when complete, provide accurate, dust- 
independent measurements of the SFRs in our redshift sam- 
ple. The ongoing SHADES survey of bright, F(850^m) > 
8m Jy, sub-mm galaxies will provide a measurement of the 
clustering of these probable progenitors of the ERGs (Dun- 
lop 2005; Mortier et al. 2005). We have planned further and 
more extensive spectroscopic surveys, including in the near- 
IR, of high redshift galaxy clusters and large-scale structure. 
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